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Abstract 

Several assemblies of guanine molecules are investigated by means of first-principle calculations. 
Such structures include stacked and hydrogen-bonded dimers, as well as vertical columns and 
planar ribbons, respectively, obtained by periodically replicating the dimers. Our results are in 
good agreement with experimental data for isolated molecules, isolated dimers, and periodic rib- 
bons. For stacked dimers and columns, the stability is affected by the relative charge distribution 
of the vr orbitals in adjacent guanine molecules, vr-vr coupling in some stacked columns induces 
dispersive energy bands, while no dispersion is identified in the planar ribbons along the connec- 
tions of hydrogen bonds. The implications for different materials comprised of guanine aggregates 
are discussed. The bandstructure of dispersive configurations may justify a contribution of band 
transport (Bloch type) in the conduction mechanism of deoxyguanosine fibres, while in DNA-like 
configurations band transport should be negligible. 

PACS numbers: PACS numbers: 73.22.Dj,71.15.Mb,87.14.Gg 
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I. INTRODUCTION. 



The twofold issue of scaling down the electronic devices and of realizing a high circuit 
integration density on a single chip has recently seen the rise of the field of molecular elec- 
tronics, which consists of the use of molecules to realize electrically conducting structures 
[jl], ||. Because of their sequence-specific recognition properties, DNA molecules are attract- 
ing attention for the construction of nanometer scale devices, where they might be used, by 
virtue of their self-assembling capabilities, to wire the electronic materials in a programmable 
way 1^. This research path has led recently to a set of controlled experiments for the direct 
measurement of the d.c. conductivity. 

Using interdigital electrodes, anisotropic conductivity was found in an aligned DNA cast 
film: At room temperature, a large ohmic current, linearly increasing with the applied 
voltage, was measured [Q. Ohmic behavior and high conductivity were found, also, for a 
a 600 nm long rope made f a few A-DNA molecules Instead, nonlinear current /volt age 
curves, exhibiting a voltage gap at low applied voltage, were measured through a single 
10 nm long poly(dG)/poly(dC) DNA molecule trapped between two metal nanoelectrodes 
0. Large currents were observed, in air and in vacuum, both at ambient temperature and 
at 4 K. The authors suggested that the observed electron transport is best explained by a 
semiconductor-like band model where the electronic states are delocalized over the entire 
length of the base pair stack 0. 

The guanine (G) base is particularly interesting in view of obtaining conductive molecular 
aggregates, because of its low ionization potential, which suggests its viability to mediate 
charge motion along a sequence of bases 0, |^. This peculiar property has opened the 
way to the measurement of the the electrical conductivity of G aggregates. Interestingly, a 
metallic nanogate filled with a dried solution of a lipophilic derivative of 2'-deoxyguanosine 
1^, |T^, displayed a current /volt age behavior |11|] similar to that of the semiconducting 
poly(G)/poly(C) sample [0. Previous investigations had shown that, in organic solvents, 
this molecule undergoes extensive self-assembly, mediated by H-bonding among the guanine 
bases, to give ribbon-like aggregates that, upon drying, form fiber structures within which 
the guanine cores of the ribbons lie on parallel planes at a distance of about 3.4 A [ll2| , |13|. 
While a direct comparison among the findings for the deoxyguanosine materials and for 
DNA cannot be made, because of the different overall experimental settings and chemical 
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nature of the molecules, one cannot avoid to notice the similar and peculiar current/voltage 
characteristics of the nanogates interconnected by molecules featuring self-assembled |P, |10| 
or inherent guanine stacks. 

In the following, we have chosen to perform ab-initio calculations of the structural and 
ground state electronic properties of extended model structures whose building block is the 
G base alone. Of course, we are well aware that these model structures are only partially 
related to the structure of a poly(G)/poly(C) duplex p, where each guanine is H-bonded 
to a cytosine in the opposite strand. However, self-assembled G structures are among the 
simplest base aggregates characterized experimentally |T2|, thus allowing accurate theo- 
retical calculations as well as comparison with experimental data. Additionally, our choice is 
motivated by the role played by guanine, the DNA base with the lowest ionization potential, 
in the mechanism of charge transport [14, 15, |T^ . 

The main question that we address here is whether the electronic properties of extended 
G-based structures can account for a band-like mechanism for charge transport. Previous 
theoretical studies for DNA bases and for their assemblies can be found in the existing 



literature JT^ |T8|, |T9], |20|, ^ Both MP2/Hartree-Fock and density-functional-based cal- 
culations were performed in the past with localized basis sets (GAUSSIAN) to treat isolated 



guanines and small clusters of guanines |17|. Here, we focus on density-functional-based cal- 
culations using a plane wave basis and ab-initio pseudopotentials ||2^ : this technique should 
allow for a correct description of solid aggregates of G's, resembling those of long-range 
deoxyguanosine fibers, provided the single molecules and the clusters are well described. 
Before simulating the structures of our interest, which are planar ribbons and stacked se- 
quences of G's, we check that our technique is well suitable to describe isolated G molecules 
as well as pairs of molecules in different (lateral or vertical) configurations. We then con- 
sider (GG)„ vertical stacks, and (GG)„ isolated and stacked periodic ribbons. We study 
the stability of different vertical configurations as a function of the relative atomic positions 
between adjacent G molecules, and identify the role of tt-tt coupling. We demonstrate that 
pseudopotential plane-wave Density Functional Theory (DFT) calculations at the current 
level of accuracy are able to reproduce not only the equilibrium length of isolated hydrogen 
bonds, but also periodic sequences of such bonds. Concerning the electronic properties, 
all the model solids that we considered are semiconducting with large energy gaps. The 
dispersion of the highest valence and lowest conduction bands is always negligible in the 
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{x,y) plane containing the molecules (the guanines are connected laterally by H bonds). 
Conversely, the dispersion along the z direction perpendicular to the G planes found to be 
extremely sensitive to the detailed geometric stacking of the bases, which in turn can be 
affected by different environments. In vertical stacking geometries that are similar to those 
present in DNA, the calculated dispersion of both the valence band maximum (VBM) and 
conduction band minimum (CBM) are extremely small. Consequently, both electrons and 
holes introduced by "doping" or photoexcitation would have very large effective masses - 
hence low mobility - in these types of structures. Instead, our results show that a band-like 
contribution may be responsible for the conduction mechanism in the 2'-deoxyguanosine 
lipophilic derivative . 

The paper is organized as follows. Section II describes the computational method. Section 
III deals with the calculated equilibrium geometries and the electronic structure for the 
isolated G molecule and for the model guanine assemblies. Section IV presents a discussion 
of our results with an outlook to the implications for the physics of guanine ribbons and 
DNA-like stackings. Finally, section V contains a summary of the arguments presented in 
the paper. 

II. METHOD. 

Our calculations are based on the DFT in the Local Density Approximation (LDA) |^ . 
For the hydrogen-bonded pairs and ribbons, we take into account BLYP gradient correc- 
tions to the exchange-correlation functional [^. The electron-ion interaction is described 
via ab-initio norm conserving pseudopotentials in the factorized form of Kleinman and By- 



lander The search for optimized metastable structures is performed by total energy 

minimizations with respect to the ionic and electronic degrees of freedom. The former are 
represented by the ionic coordinates, the latter by the electronic wavef unctions. The ions 
are treated in a classical formalism, and are displaced according to the forces derived from 



the potential determined by the full quantum mechanical electronic structure |2^, within 



a Car-Parrinello-like scheme [^]. For any selected geometry, all the atoms are allowed to 



relax, until the forces vanish within an accuracy of 0.05 eV/A. Thus, for each metastable 
structure, we obtain both the geometry and the consistent single-particle electron energies 
and wavef unctions. The electron wavefunctions are expanded in a basis of plane waves with 
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kinetic energy up to 50 Ry. This cutoff is very high with respect to standard first-principle 
calculations for most solids, and is due to the presence of the first-row elements C, N, and 
O, whose valence electron wavefunctions have strong oscillations in the region around the 
nucleus, thus needing many plane waves for an accurate treatment. For one of the model 
systems, we have verified that increasing the precision of the calculations both in the plane- 
wave kinetic energy cutoff (up to 60 Ry) and in the accuracy within which the atomic forces 
vanish (0.025 eV/A), gives changes of the bond lengths smaller than 0.01 A and of the 
bond angles smaller than 0.5 degrees, within a guanine molecule, whereas the inter-planar 
distance in the stacks is not affected at all. Additionally, for the diatomic molecules N2 and 
O2 the employed pseudopotential was tested up to a cutoff of 80 Ry, with no significant 
improvement. 

Our method employs periodic boundary conditions in three dimensions. In order to 
simulate isolated molecules, we choose a large supercell of size 15.9 A x 15.9 A x 10.6 
A. Such a choice ensures that the minimum distance (in any spatial direction) between 
two molecules is larger than 8.5 A. In particular, the distance is 10.6 A in the direction 
perpendicular to the plane of the molecules, very large with respect to the distance of 3.37 
A between two neighboring bases in B-DNA. The supercells for different assemblies of G 
molecules are described in section III, where the model structures are discussed in detail. 
For Brillouin Zone sums, the single high-symmetry F point has been employed in the case 
of isolated molecules and dimers, while one or two (depending on the symmetries) special k 



points have been employed in the case of periodic columns and ribbons, and in the case 
of stacked ribbons. 

The computational technique has been successfully applied in many investigations of the 
structure and electronic properties of inorganic and organic materials ^ . In section 

III, before reporting the results of our calculations for the model guanine assemblies, we 
address the issue of extending such calculations to biomolecules, by presenting a test on the 
G base. 



III. RESULTS. 

In this section we present our results about the isolated G molecules and about their 
assemblies and discuss them in the frame of the existing theoretical literature, which is 
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limited to isolated G's and G-pairs |jT8|, |T9|, ^ and does not give any account of periodic 



G columns and ribbons. We relate the outcome of the calculations to new experimental 
findings about the structure and electrical behavior of deoxyguanosine-based solids [1T0| . 

The section is divided into sub-sections for the different systems: (A) the isolated G 
molecule; (B) the stacked GG dimers; (C) the stacked (GG)„ columns; (D) the hydrogen- 
bonded GG pairs; (E) the hydrogen-bonded (GG)„ ribbons and the stacked ribbons. 

A. The Isolated Guanine Molecule. 

We start by presenting the calculated structure and selected electron states for the G 
molecule. This simple system allows us to evaluate the accuracy of our method, by compar- 
ison with X-ray data and with the outcome of previous ab-initio calculations. Furthermore, 
the total energy of the equilibrium structure of the isolated molecules is necessary in order 
to evaluate formation energies of dimers, ribbons, and columns. 

The structure of guanine is well known from X-ray studies |31], Figure |l](a,b) shows 



a plane view of the molecule, and the isosurface plot for the total charge density. The cal- 
culated bond lengths and angles are in good agreement with X-ray data (see Table |). Most 
of the bond lengths are underestimated within 2% and 3% with respect to the experimental 
structure, the only exceptions being the underestimate of 4.3% for the Cq-Oq bond and the 
overestimate of 1.6% for the Cg-Ni bond. The average C-H and N-H distances, not reported 
in Table I, are 1.08 A and 1.01 A, respectively, in good agreement with bond lengths in NH3 
and CH3 molecules. The bond angles in guanine are also reproduced with a high degree of 
accuracy (Table H): the discrepancy with respect to the experimental data is below 2.5%, 
but the average percentage error is 1%. The formation energy of G, with respect to its 
elemental components in stable phases (O2, H2, and N2 molecules, crystalline diamond) is 
87.5 Kcal/mole. By starting the atomic relaxation with a planar G molecule as the initial 
condition, we find that the planar configuration is indeed a metastable state. By consider- 
ing a different initial condition with a non-planar amino group (the NH2 complex bonded 
to the site C2, see Figure |l|a), we find another metastable state. However, the deviation 
from planarity is small, and the total energy difference between the planar and the puckered 
geometries is also very small, within the precision of our calculations (estimated to be about 
10 meV/G P3[). While previous calculations have pointed out a stronger stabilization effect 
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of non-planarity |T^, no direct gas phase experimental data exist. Moreover, the hydrogen 
bonds tend to flatten the structure when forming dimers and ribbons of G's; thus, we are 
confident that our results for the periodic structures are not affected by this issue. 

The single particle eigenvalues are characterized by a DFT-LDA energy gap of 4.8 eV 
between the highest occupied (HOMO) and the lowest unoccupied (LUMO) electron states 
2^ . We have identified cr and tt orbitals. The HOMO (see Figure |l|c) has a tt character and 



is localized on the Cg, Og, N2 atoms, and on the C4-C5 and N3-C2 bonds. The LUMO (see 
Figure |l|d) has a tt character and is localized on the Ng, C4, C5, C2, and Ni atoms. Because 
they extend out of the guanine plane, both the HOMO and the LUMO states are well 
suitable for interactions with adjacent similar states when forming vertical stacks, inducing 
a splitting of degenerate molecular orbitals by an amount that depends on the strength of 
such interactions. In the case of an infinite periodic stack, this is the mechanism that may 
give rise to band dispersion for sufficiently strong coupling, and to mobile carriers if the 
bands are partially occupied (for instance, as a consequence of doping or photoexcitation). 



B. Stacked GG dimers. 



In order to select low-energy geometries for the vertical columnar structures, whose elec- 
tronic properties are the main subject under investigation, we have first considered stacked 
dimers. We name stacked dimer a pair of G molecules lying in parallel planes whose distance 
in the perpendicular direction is an output of our calculations. A column is obtained by 
periodically replicating a dimer along the stacking (perpendicular) direction. 

We have analyzed several configurations (Figure |^) characterized by the relative azimuthal 
rotation angle of the two G's in the pair, with respect to an axis perpendicular to the G 
plane, and by an in-plane translation. The supercell used in the calculations was 15.9 A x 
15.9 A X 19.1 A. With periodic boundary conditions, two dimers in neighboring supercells 
are 15.7 A apart in the stacking direction, sufficient to avoid spurious interactions between 
them. By allowing all the atoms to relax, the interplanar distance between the two G's in a 
pair was a free parameter. In this way we determined the equilibrium interplanar distance 
to be used in the calculations for the periodic columns. For all the configurations shown in 
Figure (e.g., independently of the relative rotation angle between the bases), the average 
interplanar distance was 3.37 A, typical of base stacks in B-DNA. Each base maintained a 
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planar geometry, with out-of-plane fluctuations smaller than 0.05 A. 

Top views of the computed stacked dimers are illustrated in Figure |^. In Figure the 
two G's of the dimer are perfectly eclipsed, in particular the hexagonal rings are on top of 
each other, and there is maximum superposition of the vr-like HOMO and LUMO of the 
two guanines (label GGv.A). In Figure ^b, the rotation angle is zero as for GGv.A, but the 
center of mass of one molecule is shifted with respect to the other, in order to lower the 
TT-Tc superposition (label GGv.B). In Figure |^c, the azimuthal rotation angle is zero, as well 
as the relative translation of the two G's, but there is a reflection of the upper molecule 
with respect to its plane: this configuration (label GGv.C) is the only one, among those 
considered in this work, that exhibits a reflection and allows to discriminate between two 
molecular faces. In Figure ||d, the rotation angle is 180°, and the tt-tt superposition is large, 
though smaller than in GGv.A (label GGv.D). In Figure 0e, the azimuthal rotation angle is 
36°, and a translation brings the two molecules in a configuration similar to that in B-DNA 
(label GGv.E). 

Although some of these stacked structures are not likely to occur in nature, they allow 
to understand important microscopic features that may be relevant in real structures. Most 
notably, the dependence of the stability and of the electronic properties on the azimuthal 
angle and vr superposition is accessible. The lowest-energy configuration among those shown 
in Figure ^ is GGv.D, while GGv.A has the highest formation energy, and the other dimers 
have intermediate formation energies |^|. However, the energy difference between the two 
extreme cases GGv.A and GGv.D is small, about 250 meV/G. We attribute the highest 
formation energy of GGv.A to the electrostatic repulsion due to the complete tt-tt superpo- 
sition. In fact, the yr-like HOMO's of guanines in neighboring planes are mostly responsible 
for their interaction: the superposition of negative charge in the same region of space (the 
hexagonal ring) for configuration GGv.A contributes with a Coulomb repulsion. The tt-tt 
superposition is large also in GGv.D, but the repulsion is much smaller, making GGv.D a 
more viable model for a stacked GG pair . 

We wish to point out that, although the lower superposition of adjacent tt orbitals 
of GGv.D with respect to GGv.A decreases the electrostatic repulsion between the two 
molecules and makes the dimer more viable, such interaction is still rather large and is the 
origin of energy dispersion in columnar structures based on the GGv.D dimer. This issue 
is addressed in the next sub-section, by computing the band structures of model periodic 



8 



stacks built up with the dimers just described. It is also worth stressing that, for the stacked 
dimers, without inclusion of sugars, phosphates, and water as in real situations, the rotation 
angle of 36° is not preferred. This result is in line with the fact that the inner core of the 
base pair stack is not completeley responsible for the stability of double-stranded DNA, but 
also the backbone and the environment are relevant factors. 



C. Periodic stacked poly(G) columns. 

The building block of each periodic stack is a GG pair; we have analyzed five different 
model columns starting from the dimers described above (Figure ^). The supercell for these 
calculations is 15.9 A x 15.9 A x 6.74 A: no vacuum region separates two adjacent dimers 
in the stacking direction. The periodic structures are labeled with the same names as the 
stacked dimers. Direct and reciprocal one-dimensional crystal lattices are associated to the 
periodic columns: the basis vectors of these lattices are a3=6.74 A, and b3=||a3. 

The formation energies are reported in Table 0. The energetical order is the same as 
for the dimers. The most stable configuration is GGv.D. We stress again that configuration 
A, with very large tt-tt superposition, has the highest formation energy. The energy is 
substantially reduced in the column GGv.D, where the superposition is still large but the 
repulsion is weaker: the hexagonal rings lie on top of each other, but the O atoms in adjacent 
planes lie opposite to each other. Thus, we attribute the dominant repulsive contribution of 
the TT-TT interaction to the O atoms. 

A detailed analysis of the electronic properties reveals other interesting features. We 
report the numerical data in Table |T| and we show the bandstructure for columns GGv.A, 
GGv.D, and GGv.E in Figure |. 

GGv.E, with the two G's rotated by 36° as in portions of nucleic acids, has fiat HOMO- 
and LUMO-derived bands, very large effective masses, and is therefore incompatible with 
even a partial band transport mechanism for the conduction in guanine-based devices [|10|. 



The periodic columns GGv.A and GGv.D, instead, have electronic properties that may 
support electronic transport through the base stack. For GGv.A, the HOMO-band disperses 
downwards by 0.65 eV and the LUMO-band disperses upwards by 0.52 eV, between the 
center (F) and the edge (A=ib3, b3 being the reciprocal lattice vector along the stacking 
direction) of the Brillouin Zone. For GGv.D, the dispersions of the HOMO- and LUMO- 
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bands are, respectively, 0.26 eV downwards and 0.13 eV upwards. The effective masses 
of GGv.A and GGv.D, reported in Table |I[ are of the order of l-i-2 free electron masses. 
Although these values are much larger than those of conducting organic polymers [^, they 
are similar to those of inorganic materials for which band transport is demonstrated. Notable 



examples of such materials are wide-bandgap semiconductors, such as group-Ill nitrides 
As a note to support the strength of our results, we wish to point out that a test calculation 
for the model GGv.A demonstrates that the values of the band dispersions do not depend 
(within 0.02 eV) on the accuracy required for vanishing atomic forces (0.025 eV/Ainstead 
of 0.05 eV/A), and to the approximation employed for the exchange- correlation functional 
(BLYP instead of IDA). 

To illustrate how the vr-vr stack may originate channels for charge migration through a 
band transport mechanism, in Figure |^ we show an isosurface of the HOMO state at the A 
point for the columnar structure GGv.D: the interaction between the two molecules in the 
cell is evident in the superposition resulting in a delocalized orbital. 

Summarizing our study of the columnar structures, we emphasize that the orbital inter- 
action is strong. At variance, we will soon show that it is practically absent in hydrogen- 
bonded pairs. This is in agreement with the common knowledge that hydrogen bonds have 
an electrostatic (rather than covalent) nature. In the following subsections, we compare 
the behavior of planar G pairs and ribbons with that of stacked G pairs and columns. In 
particular, we find that hydrogen bonding does not give origin to dispersive electron bands, 
contrary to what we have seen for the tt stacking. 

D. Planar hydrogen-bonded GG dimers. 

We have selected one possible arrangement of hydrogen bonds, giving a structure named 
GG3 ^7\, ^ (Figure ||a). We have not considered other documented hydrogen-bonded GG 



pairs [IT^, because we focus our interest here in the guanine ribbons present in the fibers of a 



lipophilic derivative of 2'-deoxyguanosine: Such fibers [|T3| have the bond network illustrated 
in Figure |^(b). The equilibrium structure that we obtain after atomic relaxation is in good 
agreement with that of the GG3 hydrogen-bonded pair previously described via quantum 



chemistry and DFT cluster calculations [O. The two individual G molecules in the planar 



pair remain very similar to their isolated form. The N7(II)- ■ -Nl hydrogen bond has a length 
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of 2.93 A and forms an angle of 177° (2.96 A and 173° in the theoretical literature fl^). The 
N2(H)- • -06 hydrogen bond has a length of 2.85 A and forms an angle of 169° (3.27 A and 



166° in the theoretical literature [|T3])- Our calculated value for the DFT-BLYP energy gap 
23 of the GG3 dimer is 2.45 eV. The formation energy of the structure, with respect to two 
isolated guanine molecules, is -310 meV/G. This value, compared to the formation energy of 
the most stable stacked dimer GGv.D (~ -200 meV/G), is in agreement with the previous 



demonstration that hydrogen bonds are stronger than stacking interactions ||17|| . 

The structure of the GG3 dimer was calculated only as a basic building block of the fiber- 
state ribbons, whose electronic properties we are interested in. Thus, we believe that the 
underestimate of the N2(H)- • -06 distance with respect to the results of other computations 
is not a serious issue. A slight underestimate of the hydrogen-bonds calculated in the frame 



of DFT is documented Furthermore, the contraction of the N2(H)- • -06 hydrogen bridge 
is due to an imprecise account of secondary interactions such as C8(H)- • -06. It is likely 
that such secondary interaction is reduced in one-dimensional ribbons, thus minimizing the 
shortcomings of DFT. 

E. Planar hydrogen-bonded ribbons. 

In a periodic ribbon obtained by piling up replicas of the GG3 dimer (Figure ^d), the 
equilibrium structure maintains the bonding characteristics of the dimer. There is a huge 
energy gain in forming the one-dimensional ribbon, of 820 meV/G with respect to isolated 
G molecules, and of 510 meV/G with respect to the hydrogen-bonded GG3 pair (Table p| ). 
Note that G has a large dipole moment [|19|; the dipole moments of G's add up to give a non- 
vanishing dipole also in the GG3 dimer and in the ribbon, parallel to the ribbon axis p!0| . 
Such electrostatic interactions account for the high stability of the planar hydrogen-bonded 
structures, well known both in solution and in the solid state. 

The bandstructure of the ribbon was calculated along a symmetry line parallel to its axis. 
The HOMO- and LUMO-derived bands are separated by a DFT-BLYP energy gap |2^ of 



3.84 eV, and are both dispersionless. Consequently, these bands have practically infinite 
effective masses, so that the electrons and holes in these states are not mobile according 
to this picture. The electronic state analysis shows that the HOMO and the LUMO have 
a TT character, similar to isolated G. The electronic states are localized around single G 
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molecules, no delocalized intermolecular states, extended through the hydrogen bonds in 
the ribbon, are present. The dispersion induced by hydrogen-bonding is not compatible 
with band transport. The only possible conduction mechanism through hydrogen bonds 
is via electrostatic interactions. Therefore, in a device where dried deoxyguanosine fibers 



are deposited between two metal electrodes |]T0|, if the ribbons are stretched between the 
electrodes, band transport cannot contribute to conduction. 

In order to investigate the competition of tt-tt coupling versus hydrogen-bonding interac- 
tions, we have simulated two different configurations of stacked ribbons, periodic along the 
stacking direction. Top views are shown in Figure ^(b,c). Figure ^ shows indeed a single 
ribbon, but the top view is equivalent for two exactly eclipsed ribbons on top of each other 
(label SR. A): the stacking between adjacent bases is similar to that of column GGv.A, and 
all the individual bases are perfectly aligned. Figure ^ shows a stack of two ribbons where 
only half of the individual G molecules lie on top of each other, with a stacking similar to 
that of configuration GGv.D (this structure is labeled SR.D). In both the SR.A and SR.D 
geometries, the projections of the ribbon axes on the (x, y) plane coincide, identifying an 
axis for the stack: this axis defines the T — X direction in the one dimensional BZ. 

The supercells are 21.2 A x 11.3 A x 6.74 A for SR.A and 24.3 A x 11.3 A x 6.74 
A for SR.D. The choice of different supercells allows to have a vacuum region of the same 
volume between neighboring ribbons in the two configurations SR.A and SR.D. The relaxed 
structures maintain the geometry of the single ribbon: no G deformations, no variations of H- 
bond lengths and angles, and no out-of-plane buckling are observed. This finding indicates 
that the stacking does not affect the hydrogen-bonding mechanism that determines the 
structure of the isolated ribbon. 

Numerical data for the energetics and the electronic properties are reported in Table 



ri|. Both structures SR.A and SR.D are energetically favorable with respect to isolated 



G molecules and to stacked columns, but unfavorable with respect to isolated ribbons, as 
indicated by the smaller energy gain of the stacked ribbons. The configuration SR.D is more 
stable than SR.A by 280 meV: this trend is consistent with what found previously for the 
columnar stacks, where a higher stability is achieved by lowering the tt-tt superposition in 
building the stacks. Although the configurations SR.A and SR.D that we have calculated 
here are not viable models for piling up ribbons into solid-state crystals, ordered phases for 



deoxyguanosine fibers are known to exist |12|. We note that other factors, such as presence 
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of the phosphate-sugar backbone, or the possibihty of different relative positions of the 
neighboring ribbons, should be taken into account to achieve a full description, which is 
beyond the scope of this work. 

In Figure]^ we show the calculated bandstructure for the SR.A and SR.D configurations. 
In both cases, the HOMO- and LUMO-derived bands are dispersionless along the T — X 
direction, and the charge carriers are not mobile in this direction. This behavior is the same 
as in isolated ribbons: therefore, it is a further evidence that the base-base interactions due 
to the stacking do not change the features of the ribbons. In the F — ^4 direction, parallel 
to the stacking direction, the HOMO- and LUMO-derived bands are dispersive, as in the 
columnar structures of subsection IIIC. 

The bands for structure SR.A in the T — A direction are more dispersive than those for 
structure SR.D, due to the higher tt-tt superposition. As a consequence, the corresponding 
electron and hole effective masses along the T — A direction are smaller. The effective masses 
of the stacked-ribbon periodic structures are very similar to those of the stacked periodic 
columns (subsection IIIC): this is an indication that the hydrogen-bonding network of the 
individual ribbons does not affect the electronic properties in the perpendicular direction. 
The only exception to this trend is me, which becomes infinite for the stacked ribbons SR.D 
while it is finite for the equivalent stacked column GGv.D: this is due to the fact that in 
structure SR.D one half of the individual bases are on top of each other (see Figure |^), 
while the other half do not participate in the tt-tt coupling. 

Summarizing the results discussed above, we wish to point out that the tt-tt coupling and 
H-bonding interaction mechanisms are not in competition. H-bonding has an electrostatic 
nature and accounts for in-plane stabilization: it is not compatible with a band transport 
mechanism, and it is not modified by base stacking, tt-tt coupling is weaker than H-bonding: 
it is compatible with band transport along the stacking direction and it is not affected by 
the local details of the planar base sequence. 
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IV. DISCUSSION. 



A. Guanosine films. 

As outlined in the Introduction, the recent hvely research activities in molecular electron- 
ics have pointed out peculiar guanine assemblies to be exploited as electrical conductors. 



In fact, the experiments on a lipophilic derivative of 2'-deoxyguanosine [|T0|, demonstrated 
conduction through the biomolecular material deposited in a nanogate between metal elec- 
trodes. The details of the conductivity depend on the on the experimental conditions, in 
particular on the gate width. Our results indicate that the observed conductivity, resem- 
bling that of a semiconductor in the intermediate gate length regime (few hundreds nm) and 
that of a diode junction in the short gate length regime (less than 100 nm), is compatible 
with a Bloch contribution to the transport mechanism. We have shown that delocalized 
orbitals through the base stack may be formed, provided that a relevant superposition of 
the hexagonal rings of the single guanines is maintained. Therefore, if the guanosine ribbons 
in the device gate align locally with their plane perpendicular to the direction connecting 
the electrodes, in a geometry similar to that proposed as SR.A, then extended Bloch-like 
states may be partially responsible for the observed conduction. This stacking orientation 
does not need to be complete in order for the proposed mechanism to work: it is sufficient 
that randomly aligned stacks form locally, and that the total resulting component in the 
direction connecting the electrodes is non vanishing. In the gate length regime ranging from 
100 nm to 300 nm, semiconducting-like conductivity is revealed: we propose that, in such 
a condition, the ribbons are locally stacked in such a way to form partially delocalized or- 
bitals that give a global band-like contribution. In the gate length regime below a 100 nm, 
a diode-like behavior is observed: this characteristic is assigned to an interaction of the vr 
stack with the total dipole moment of the ribbons, as discussed elsewhere |I0| . 

B. DNA molecules. 

While our study is limited to a single type of nucleoside, it also allows a discussion 
of the electronic properties of DNA molecules, where G sequences in base stacks play an 



important role |15, 3/, 38|. Depending on the energetics of the base sequence, and on the 



overall structural aspects of the system under investigation, the mechanisms proposed for 
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DNA-mediated charge migration include single-step superexcliange |3^, 0, multistep hole 
hopping P, |14|, phonon-assisted polaron hopping and band transport [§. 

In DNA double strands, the relative arrangement of neighboring bases along the axis 
of the helix is characterized by a rotation angle around 36 degrees. Our results show 
that periodic G columns in such a configuration do not support the formation of extended 
molecular orbitals. Therefore, band transport would not be effective. For what concerns 
the mechanism of conductivity in DNA molecules, our results support the conclusion that 
the contribution of band transport would be very small |]T^ ^ in contrast to recently 
proposed interpretations of experimental data , unless structural distortions, possibly ac- 
tivated by temperature effects, may induce rotations that support the formation of partially 
delocalized electron states. Although in real DNA molecules the backbone phosphate and 
sugar groups may affect the overall charge mobility, we believe that the presence of the outer 
mantle would not change our present conclusion that band transport is not supported by 
the native DNA stack: eventually, it may contribute a hopping or ionic mechanism, but it 
is unlikely to contribute to the formation of extended electron orbitals. The importance of 
the base stack for electron transfer through DNA molecules has been highlighted in recent 



theoretical invesigations O]. 



V. CONCLUSIONS. 



We have reported the results of ab-initio calculations for the structure, energetics, and 
electronic properties of several guanine assemblies. The geometry of isolated molecules and 
hydrogen-bonded dimers is well described with our technique, based on plane-wave pseu- 
dopotential density functional theory, as demonstrated by comparison to available theoretical 
and experimental data. 

We found that hydrogen-bonding and tt-tt coupling are independent mechanisms that 
control the self-assemblying of guanine bases. Hydrogen-bonding is not responsible for band 
transport: In fact, we have shown that no band dispersion is present along a planar ribbon 
of H-bonded guanines. Instead, base stacking is accompanied by tt-tt interactions that, for 
the case of sufficiently large overlap between adjacent tt orbitals {e.g., GGv.D configuration), 
induce energy dispersion and are consistent with charge mobility. Therefore, band trans- 
port may be partially responsible for charge mobility in nucleotide aggregates, in structures 
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characterized by a large base-base superposition. This mechanism is hkely complemented by 
hopping to connect (through space) different regions where such superposition is realized. 
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(b) Isosurface plot of the total charge density from the ab-initio calculation, (c) Isosurface plot of 
the HOMO, (d) Isosurface plot of the LUMO. 

FIG. 2: Top views of different stacked dimer configurations. The structures are explicitly defined 
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TABLE I: Structural data for the isolated guanine molecule. The experimental structure taken as 



reference |31, 32| is that for guanosine, with a sugar moiety instead of an H atom in position 9 
{e.g., attached to atom Ng in Figure |l|). 







distance 


(A) 


angle (°) 


this work 


experiment 


discrepancy (%) 


this work 


experiment 


discrepancy (%) 


(Cg-Nr) 


1.281 


1.304 




-1.8 


(NgCgNr) 


112.4 


113.5 


-1.0 


(C5-N7) 


1.347 


1.389 




-3.0 


(C8N7C5) 


105.5 


104.2 


+1.2 


(Cs-Ce) 


1.419 


1.415 




+0.3 


(N7C5C4) 


110.6 


110.8 


-0.2 


(Ce-Oe) 


1.186 


1.239 




-4.3 


(N7C5C6) 


130.8 


130.1 


+0.5 


(Ce-Ni) 


1.415 


1.393 




+1.6 


(C4C5C6) 


118.6 


119.1 


-0.4 


(C2-N1) 


1.345 


1.375 




-2.2 


(CsCeOe) 


131.4 


128.3 


+2.4 


(C2-N2) 


1.336 


1.341 




-0.4 


(CsCgNi) 


109.2 


111.7 


-2.2 


(C2-N3) 


1.287 


1.327 




-3.0 


(OeCeNi) 


119.4 


120.0 


-0.5 


(C4-N3) 


1.323 


1.355 




-2.4 


(C6N1C2) 


127.1 


124.9 


+1.8 


(C4-C5) 


1.383 


1.377 




+0.5 


(N1C2N3) 


122.6 


124.0 


-1.1 


(C4-N9) 


1.346 


1.377 




-2.3 


(N1C2N2) 


117.0 


116.3 


+0.6 


(Cg-Ng) 


1.361 


1.274 




-1.0 


(N2C2N3) 


120.3 


119.7 


+0.5 












(C2N3C4) 


113.9 


111.8 


+1.9 












(N3C4C5) 


128.5 


128.4 


+0.1 












(N3C4N9) 


126.0 


126.0 


0.0 












(C5C4Ng) 


104.6 


105.6 


-0.9 












(C4N9C8) 


107.0 


106.0 


+0.9 



FIG. 3: Bandstructure for configurations GGv.A, GGv.D, and GGv.E, calculated along the sym- 
metry line parallel to the stacking direction. Larger dots indicate HOMO and LUMO states. The 
single-particle energies reported in these plots are relative to the top of the highest valence band. 



FIG. 4: Isosurface plot of the HOMO state at the point A for the columnar structure GGv.D. 
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TABLE II: Energetical and electronic data of different periodic stacked configurations. Ej^^^ 
is calculated with respect to tfie N2, O2, H2 molecules, and C in the diamond phase. Ej^^^ 
is calculated with respect to isolated guanine molecules. Egap is the energy difference between 
the HOMO and LUMO single particle eigenstates, calculated by DFT-LDA [H). mg (m/^) is the 
electron (hole) effective mass, in units of the free electron mass niQ. The structures are labeled as 
in Figure ^; G is a single molecule. 





G 


GGv.A 


GGv.B 


GGv.C 


GGv.D 


GGv.E 


^orm (meV/G) 
Egap (eV) 
me (mo) 
m/i (mo) 


-3.80 


4.8 


-3.69 
-hlOO 
2.97 
1.41 
1.04 


-3.79 


3.71 


-3.93 
-130 
4.12 


-4.07 
-280 
3.63 
2.80 
2.20 


-3.89 
-100 
3.54 
00 
5.25 



TABLE III: Energetical and electronic data of isolated and stacked ribbons, as defined in Table ||. 





ribbon 


SR.A 


SR.D 


^}orm (meV/G) 


-820 


-310 


-580 


Egap (eV) 


3.84 


2.61 


3.12 


nie (mo) 




1.24 


00 


nife (mo) 




1.07 


2.05 



FIG. 5: (a) Top view of the hydrogen-bonded GG dimer. (b) Top view of the hydrogen-bonded 
ribbon obtained by periodically replicating the dimer, representing also the structure SR.A for the 
stacked ribbons, (c) Top view for the stacked ribbon structure SR.D. Black and gray dots identify 
atoms lying in two different parallel planes. 

FIG. 6: Bandstructure for configurations SR.A and SR.D, calculated along the symmetry lines 
parallel to the ribbon direction (F — X) and to the stacking direction (T — A). Large dots indicate 
HOMO and LUMO states. 
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